Spectroscopy of a fractional Josephson vortex molecule 
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In long Josephson junctions with multiple discontinuities of the Josephson phase, fractional vortex 
molecules are spontaneously formed. At each discontinuity point a fractional Josephson vortex 
carrying a magnetic flux |$| < $0, &o ~ 2.07 x 10~ Wb being the magnetic flux quantum, is 
pinned. Each vortex has an oscillatory eigenmode with a frequency that depends on c& / <&o and lies 
inside the plasma gap. We experimentally investigate the dependence of the eigenfrequencies of a 
two- vortex molecule on the distance between the vortices, on their topological charge p = 2n&/&o 
and on the bias current 7 applied to the Josephson junction. We find that with decreasing distance 
between vortices, a splitting of the eigenfrequencies occurs, that corresponds to the emergence of 
collective oscillatory modes of both vortices. We use a resonant microwave spectroscopy technique 
and find good agreement between experimental results and theoretical predictions. 

PACS numbers: 74.50.+r, 85.25.Cp 
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I. INTRODUCTION 

During the last decade the physics of fractional Joseph- 
son vortices has attracted a lot of attention^— In con- 
trast to the well known fluxons^ fractional vortices carry 
the magnetic flux of |$| < $0. Here we focus on 
the fractional vortices that emerge at the discontinu- 
ities of the Josephson phased— in long Josephson junc- 
tions (LJJ) and disregard all other types, e.g., splintered 
vortices J^r— The phase discontinuity k can be created by 
two tiny (compared to the Josephson penetration depth 
Aj) current injectors^ ' 12 ' 16 so that k oc I lni , where I lni is 
the current sent through the injectors. 

A well known property of a long Josephson junction is 
the presence of a gap in the spectrum of electromagnetic 
waves (plasma waves) — the so-called plasma gap, which 
ranges from to the plasma frequency oj p . The value of 
w p Vie is ~ 10 . . . 1000 GHz, where j c is the critical 
current density of the junction. In a LJJ with a fractional 
vortex pinned at the phase discontinuity, the vortex has 
an eigenfrequency u>q which is situated within the plasma 
gap (localized mode with a wave number k = 0). The 
value of wo(p, 7) depends on the topological charge p of 
the vortex and the normalized bias current 7 applied to 
the LJJ.il 

The topological charge of a single fractional vortex in 
an infinite LJJ with a k discontinuity of the phase is 
defined as p = 0(+oo) — <p{— 00) — k, where 4>{x) is the 
Josephson phase, that is ^-discontinuous at one point. 
In simple words, p is a continuous advance of the phase 
when one goes from x = — 00 to s = +00. Among all 
possible combinations of the phases at x = ±00 only two 
irreducible combinations are stable^: the direct vortex 
with p = — k and a complementary vortex with p = 
— k + 27rsgn(K). We assume that |k| < 2ir, otherwise the 



structure cither cannot be constructed (not a solution 
of the sine-Gordon equation) or is unstable — it splits 
into one or several fluxons and a fractional vortex with 
|p| < 2-7T.— For p — > there is a smooth transition to the 
flat phase state <p(x) = of the LJJ, while for \p\ — > 2ir 
a(n) (anti)fluxon is formed. In contrast to the fractional 
vortex this fluxon is not pinned anymore and can move 
freely.— 

The eigenfrequency of a single fractional vortex 
as a function of p and 7 was already measured 
experimentally^*^ using microwave spectroscopy. Good 
agreement with numerical calculations was found. 

Nowadays more complex fractional vortex systems, e.g. 
fractional vortex molecules consisting of two, three, etc. 
vortices can be created and studied experimentally^ In 
particular, two-vortex molecules with degenerate ground 
states^ are proposed as a candidate for a macroscopic 
qubit^S In such a molecule, the coupling between vortices 
results in shifting/splitting of their eigenfrequency^ very 
similar to the splitting of electronic energy levels in a real 
molecule when it is formed from atoms. 

In this paper we experimentally investigate the 
eigenfrequencies in two types of symmetric two-vortex 
molecules: one with parallel (ferromagnetic) configura- 
tion and the topological charges (p, p) and one with an- 
tiparallcl (antiferromagnctic) order and the topological 
charges (p, — p). We use a spectroscopic measurement 
technique^ to determine the eigenmodes of the molecule 
as a function of applied bias current and topological 
charges. 

The paper is organized as follows. In Sec.|TI]we present 
the underlaying theory and present numerically calcu- 
lated ground states of molecules as well as their eigen- 
frequencies. Experimental results are presented and dis- 
cussed in Sec. IHII Section IIVI concludes the work. 
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II. THEORY 

The perturbed sine-Gordon equation, which describes 
the dynamics of the Josephson phase is given by-i£ 

4>xx ~ <f>tt - sin(» = a<j) t - 7(2;) + 9 xx (x), (1) 

where subscripts x and t denote derivatives with respect 
to space and time and the step function 0(x) describes the 
position of discontinuities. In the case of one re disconti- 
nuity located at x = 0, 9{x) = kH(i), where H(x) is a 
Heaviside step function. In Eq. ([T]) the coordinate x and 
time t are normalized to the Josephson penetration depth 
A j and the inverse plasma frequency w" 1 , respectively. 
a = \j\fJTc is the dimcnsionlcss damping parameter, j3 c 
is the Stewart-McCumbcr parameter. The bias current 
density 7 is normalized to the critical current density j c 
of the junction. To solve Eq. (TTJ) numerically we intro- 
duce a new continuous phased n(x,t) = 4>{x,t) — 6(x) 
and solve the resulting smooth equation for fx(x,t). 

The eigenfrequency of a single vortex with a topologi- 
cal charge p in an unbiased (7 = 0) infinite LJJ without 
damping (a = 0) is given byi£ 



w (p) =w p W-cos| fcos^ 



»s 2 K\. (2) 



In case of 7 7^ no analytical expression has been found 
and therefore the eigenfrequency has to be calculated nu- 
merically. A good approximation to the numerical results 
is given by 



from each other, where I = 2ttR/Xj is the mean cir- 
cumference of the junction normalized to Xj. There are 
certain advantages of this particular geometry. First, due 
to a particular shape of the electrodes, annular junctions 
have intrinsically a much more homogeneous bias cur- 
rent density than linear LJJs^ which is important to 
measure uo(j). Second, any real LJJ has a finite length, 
which results in a somewhat modified phase profile and 
eigenfrequency spectrum of the vortex molecules. The 
boundary conditions of the junction thereby impose an 
additional level of symmetry /ordering on the molecule. 
In an annulus this additional symmetry is always the 
same as the internal configuration of the molecule, FM 
or AFM^ 

In an annular junction each vortex has the same neigh- 
bor on both sides (right and left) so that magnetic field 
profiles and eigenfrequencies are somewhat different from 
the results obtained earlier for an infinite LJjJ^ This is 
especially true for short JJs. In principle, there are many 
eigenmodes ujq < uj\ < . . . enumerated according to their 
frequency. In the following we use w+ to denote the eigen- 
frequency of the vortex in-phase oscillations, while oj_ de- 
notes the oscillation mode with a phase shift of w. They 
correspond to the two lowest eigenmodes loq and a>i, i.e. 
lu± = ujq or Depending on the configuration of the 
molecule w+ can be larger or smaller than w_. In gen- 
eral, the eigenfrequencies of a molecule can only be calcu- 
lated numerically^ Due to the vast number of possible 
molecule configurations we, in the following, limit our- 
selves to the symmetric configurations and discuss them 
in more detail. 



wo(p,7) - w (p, 0) 
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|p| < 2^ 



(3) 
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is the depinning curren t 21 ' 22 of the vortex. For bias cur- 
rents exceeding 7 C , the fractional vortex flips into its com- 
plementary partner, thereby emitting a fluxon. 

If one considers two discontinuities with the same ab- 
solute value of |re| that are situated at x± t 2 = ±a/2 (a 
is the distance between the discontinuities) in an infinite 
LJJ, four irreducible vortex configurations are possible^ 
In the case of equal discontinuities, i.e. (—re, —re), there 
exists the symmetric parallel (ferromagnetically) ordered 
(FM) vortex molecule with (pi, p 2 ) = (+re, +re) and the 
asymmetric antiparallel (antiferromagnetically) ordered 
(AFM) vortex molecule with (pi,p2) = (+re, +re — 27r). 
For discontinuities of different signs, i.e. (+re, — re), 
one gets either the symmetric antiparallel molecule with 
(pi,p2) = (—«,+«) or the asymmetric parallel vortex 
molecule with (pi, P2) = (— re + 27r, +re). 

In the experiments described below we use annular 
LJJs 2 ^ with two discontinuities at a distance a = 1/2 



A. Parallel (FM) molecule 

The magnetic field oc fi x (x) of a molecule consisting of 
two parallel ordered vortices pinned at (—re, — re) discon- 
tinuities is exemplarily shown in Fig. [Ha) f° r an annu- 
lar LJJ of length 1 = 6. The vortex distance is half of 
the junction length. No bias current is applied. Then, 
the FM configuration is stable for |re| < 2ir and for all 
lengths due to annular boundary conditions and a = 1/2, 
cf. Refs. Il7l and [l8l for more general cases. For \k\ > 2n 7 
a single direct vortex with |p| > 2-7T is not stable any- 
more. Instead, at each discontinuity a(n) (anti)fluxon 
and a (re — 27rsgn(re))-vortex (all of the same polarity)^ 
is formed. 

The eigenfrequencies ui± of a FM ordered molecule are 
shown in Fig. [21 a) as a function of re for different vortex 
distances. Here, uj + < w_. The dependence of 0J±(a) is 
shown in Fig. [3] for two molecules with specific values of 
the discontinuities re^ = re. For comparison, the eigenfre- 
quency of a single vortex in an infinitely long JJ is also 
shown. As can be seen the frequency splitting increases 
for smaller a (larger coupling). However, the divergence 
of the out-of-phase mode u;_ = u\ for a — > is exclusively 
related to the fact, that also I = 2a — s- and therefore 
lui cx l/l — > 00 for I < 1. 
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FIG. 1: Magnetic field profile fj, x (x) in an annular LJJ of 
length I — 6 and a — 3 for different values of k: (a) parallel 
configuration, (b) antiparallel configuration. 



B. Antiparallel (AFM) molecule 

Figure QJb) shows the magnetic field oc pL x {x) for the 
antiparallel molecule pinned at (— ft, +k) discontinuities 
in an annular LJJ of I = 6. In contrast to the parallel 
configuration there exists a critical value k c (tt < k c < 
27r), above which the vortices are not stable anymore: 
after exchanging one 4> the (+ft, — ft)-molecule turns into 
a complementary one (+ft — 2ir, —k + 2tt). The value of 
k c depends on a and can in general only be calculated 
numerically. For an AFM molecule in an infinitely long 
J J k c was calculated in Ref.[l8l Note that for (2-7T — k c ) < 
\k\ < n c both states (+ft, —ft) and — 2tt, — k + 27t) are 
stable and that for k ^ tt these solutions have different 
energies. 

The eigenfrequencies w±(k) of the two lowest modes 
of the antiparallel configuration are shown in Fig. HJb) . 
The vertical lines describe the critical discontinuity k c for 
different vortex distances a. In contrast to the parallel 
configuration, here lj + > w_ . 

The dependence of the eigenfrequencies of two antipar- 
allel vortices for fixed k values as a function of the dis- 
tance a is shown in Fig. [3] As with the parallel configu- 
ration the frequency splitting increases when a decreases. 
One has to recognize, that for ft = 1.27T the critical vor- 
tex distance is a c rj 2.15. Therefore, the vortex state 
(+ft, —ft) becomes unstable at a — > a c and the frequency 
w_ = luq — > 0. For a < a c , uj± correspond to the comple- 
mentary state (k — 2-7T, 2ir — ft), which is stable for all a. 
In comparison to the lowest eigenfrcquency of the paral- 
lel configuration, this mode has a smaller frequency and 
this difference increases as a — > 0. 



FIG. 2: (Color online) Numerically calculated lowest eigen- 
mode(s) of a fractional vortex molecule in an annular LJJ for 
different distances a between the vortices: (a) parallel config- 
uration, (b) antiparallel configuration. The junction length I 
is twice the distance a. 

For a biased molecule (7 7^ 0) the approximation J3]) 
can also be used to estimate u>±(p, 7). 



III. EXPERIMENTS 

For experiments we used underdamped long Nb-AlOz- 
Nb annular tunnel junctions fabricated using standard 
trilayer technology^ Each JJ was equipped with two 
pairs of current injectors as shown in Fig. |4j Note that 
injector pairs are placed opposite to each other and di- 
rectly on top/below the bias electrodes of the junction, 
as earlier measurements showed this layout to give best 
agreement of the critical current dependence J c (/i n j) with 
theory, i.e. Eq. As the JJs have been fabricated in 
an overlap geometry, their static and dynamic proper- 
ties are modified by the overlap region of top and bot- 
tom electrode, the so-called idle regioni^Zr— Idle region 
effects can be taken into account by a renormalization 
of the junction length, i.e. by introducing an effective 
Joscphson penetration depth A,/. c ff > Xj. However, we 
find that the usual correction^ overestimates the idle 
region effect. Table Q] summarizes the properties of the 
measured junctions, listing Joscphson penetration depths 
and normalized junction lengths, with and without idle 
region corrections, as well as the normalized lengths we 
found to give best agreement with theory. 

The measurements shown in this paper have been per- 
formed at T — 4.2 K with the resonant vortex excitation 
technique already used to determine the eigenfrcquency 
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FIG. 3: (Color online) The two lowest eigenfrequencies oj±(k) 
of a symmetric vortex molecule in an annular LJJ with length 
I — 2a for different distances a between the vortices: (a) \n\ = 
0.6n and (b) \k\ — 1.2n. Solid (dashed) lines correspond to the 
FM (AFM) configuration. For comparison the eigenfrequency 
of a single vortex in an infinitely long JJ is also shown (dotted 
lines). Full (open) symbols correspond to fitted values of ojo 
(oil) from experimental data, see text. 




FIG. 4: Optical image (top view) of sample 2: annular LJJ 
with two pairs of injectors. The dc bias I is superimposed 
with a microwave induced ac current 7 ac , that is capacitively 
coupled into the bias leads. Imj.i is applied to the left (top) 
pair of injectors, 7mj,2 to the right (bottom) pair respectively. 



of a single fractional vortexii In brief, we continuously 
apply microwaves with fixed frequency u; oxt < uio(p) 
and fixed power to the junction while ramping up the 
bias current 7. As u>i(p, 7) decreases with increasing 7, 
cf. Eq. ([3]), the resonance condition wo(p, 7r) — w cx t is 
achieved for some j r . Then, depending on applied mi- 
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TABLE I: Sample properties measured at T = 4.2 K. Junc- 
tion widths w and injector distances Ax are 5 /an, injector 
widths Aw are 5 /jm, except for sample 3, where Aw = 2 fj,m. 
l e g = 27T_R/Aj, e ff and h = 2nR/Xj are the normalized junc- 
tion length with and without idle region corrections and R 
is the radius of the junction. I is the length which gives best 
agreement between measurement and theory. uj p is the plasma 
frequency obtained from spectroscopic measurements. 



crowave power and the internal damping of the junction, 
the vortex molecule may flip and switch the junction into 
its resistive stated If, however, the junction remains in 
the zero-voltage state it switches to the non-zero volt- 
age state by means of thermal escape when the bias ap- 
proaches the fluctuation- free critical current j c . By re- 
peatedly measuring the switching current and system- 
atically scanning microwave power and frequency we can 
map out the resonances corresponding to the eigcnmodes 
of the molecule. As the J J itself is a non-linear oscillator, 
resonances shift towards lower frequencies with increas- 
ing microwave power.— i22rJ2£ This is especially relevant at 
small 7 where the energy barrier between the zero and the 
non-zero voltage state is still large and high microwave 
power is necessary for switching. We therefore determine 
7 r from the measurements with the lowest power to keep 
non-linear effects to a minimum. 

Before spectroscopy was performed, injectors were cal- 
ibrated. For this we first measured the critical current I c 
of the junction as a function of the current through each 
injector pair ii n j,i and ii n j,2 individually. From Eq. (UJ 
and k oc linj we then determined the injector currents 
7^° necessary to create a discontinuity of Ki = 2n, where 
Ic(Iinj.i) have the cusp like first minima. Figure Eta) and 
(b) show the measured I c (Iinj,i) dependences of sample 
1, as well as the fitted theoretical curves. Although both 
injector pairs have nominally the same width and dis- 



tance, the /^j™ are usually somewhat different. In sample 
3 the inequality was introduced intentionally, as we de- 
liberately cut the top wiring layer between the injectors 
to increase the coupling of /i n j,i to the junction. How- 
ever, in samples 1 and 2, both injector pairs are identi- 
cal by design. Here, we believe, the deviations are due 
to lithographical inaccuracy and different thicknesses of 
the top and bottom Nb layer. To create symmetric vor- 
tex molecules we henceforth applied the injector currents 
Iinjj = ±K|/;™|/27r simultaneously. 

Figure [5jc) shows the measured and the numerically 
calculated 7 C = I c (k)/I c {0) dependence of a symmetric 
FM vortex molecule pinned at (k, k) discontinuities. As 
can be seen there is a good agreement between experi- 
ment and theory. Note that j c (+k, +k) looks qualita- 
tively similar to 7 C («0 for a single vortex^ yet the first 
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discontinuity k [ti] discontinuity |k| [tc] 

FIG. 5: (Color online) Measured (symbols) and numerically 
calculated (red line) critical current y c — I c (li n j,i)/Ic(0) as a 
function of injector current. Single vortex cases with either 
(a) injector current 7i n j,i or (b) 7i n j,2 applied, (c) FM and (d) 
AFM molecule cases corresponding to (+ft, +k) and (+«;, — «) 
discontinuities. Data is shown for sample 1. For the vortex 
molecule injector currents Imj.i = ±K|/ i ™j°|/27r were applied 
simultaneously. 



minimum positions are slightly different: 7™ m of a sym- 
metric FM vortex molecule in an annular JJ with a finite 
width Ax, Aw of injector pairs corresponds to a Ki value, 
which is slightly larger than 2tt. 

The 7c(+k, — k) curve for the AFM molecule is de- 
picted in Fig. [5jd). Again, experimental data and nu- 
merical calculations are in good agreement. In con- 
trast to the FM configuration, here the 7 c -minima cor- 
respond to |«| ~ (2n + 1)tt instead of k ~ 2nn, cf. 
Eq. (QJ. Note that the 7 c (+k, — k) curve is identical 
to the 7 c (k) dependence of a linear LJJ of length 1/2 
with a single central vortex^ As already mentioned in 
Sec. Ill BL the AFM molecule exhibits regions of bista- 
bility of the direct and the complementary vortex states 
around \n\ = (2n + l)ir. Between 2tt — K c (a) < k < n c (a) 
the critical current of both vortex states can be traced 
out^ see Fig.[5][d). In AFM configuration the j c (+k, —k) 
pattern sensitively depends on junction length as well as 
on injector size^i: For infinitesimally small injectors, the 
7 c (+k, — k) is 27r periodic with maxima at k = 2im and 
minima at n — (2n + l)7r. For finite size injectors, how- 
ever, this is not the case: 7 c (±27m, ^f2nn) < 7 C (0, 0) and 
the minimum positions shift to k > (2ra + l)7r. In the FM 
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FIG. 6: (Color online) Comparison of measured resonance 
peak positions (symbols) with the numerically calculated 
cjo(k, 7) dependences (solid lines). Data is shown for (a) the 
plasma frequency scan (k = 0), (b) the FM (— k, —k) and (c) 
the AFM (— k, +k) molecule configurations in sample 1. Sym- 
bols of same color correspond to the same k value. Symbols 
are explained in the text. 



case basically only I^f\ depends on injector size 



the 



shape of the curve does not change. By carefully compar- 
ing the 7 c (/i n j j i)// c (0) curves of the AFM configuration 
with numerical simulations we obtained the normalized 
junction length I, sec Tab. HI which we used for the sub- 
sequent calculation of eigcnfrcqucncies. 

After calibration of injectors we spectroscopically 
mapped out the well-known dependence of the plasma 
frequency 0^0(0,7) = uj p \/l — 7 2 on applied bias cur- 
rent. Figure IH^a) shows excmplarily the results of the 
plasma frequency scan of sample 1. All measured res- 
onance peak positions are plotted. Beside the desired 
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plasma resonance (squares), additional resonances were 
observed. They do not correspond to a sub- or a su- 
pcrharmonic excitation of the plasma resonance, but can 
be subsumized into two groups. One set of resonances 
(triangles) occurs at high frequencies and close to the 
critical current "f c of the junction. These resonances are 
most likely the result of a microwave enhanced thermal 
escape. The second group of resonances (circles) seem 
to be erratically distributed over the full bias current 
range. Yet, a comparison with other spectroscopic mea- 
surements, e.g. see Fig. &b), shows that they appear 
in rather small frequency bands, where the coupling of 
the microwave current to the junction is exceptionally 
high. We therefore believe that they are the result of a 
parasitic coupling of the junction to cavity modes of the 
sample box and the bias circuitry. Note, that in contrast 
to a typical spectroscopic measurement of the plasma fre- 
quency or the eigenfrequency of a single fractional vortex, 
e.g. see Ref. 0, we did not simply trace out the desired 
resonance curve by optimizing the measured frequency 
range and the applied microwave power. Instead, we sys- 
tematically scanned over a large microwave power and 
frequency range to ensure, that no resonance is missed. 
Evidently, this measurement technique also reveals par- 
asitic resonances, which are not of interest and must be 
removed in a manual post-selection process. From the 
measured plasma resonance (squares) we extracted to p 
and the fluctuation-free critical current Iq = j c 2-KRw of 
the junction, which then served as a reference for the 
following measurements on vortex molecules. R is the 
junction radius, w the junction width. Wc also deter- 
mined Up from measurements of the zero field steps^ on 
the current-voltage characteristics of the junction. After 
corrections due to idle region effect o 28 ' 30 ' 31 we typically 
find both values in agreement by better than 90%. 

Figure [HJb) shows all measured resonance peak posi- 
tions of a set of FM ordered vortex molecules in sam- 
ple 1. Also plotted are the corresponding numerically 
calculated u)q(k, 7) dependences. As sample 1 is rather 
short (I = 5.1), u>i = uj- is much larger than the high- 
est frequency of our rf-generator and could therefore not 
be measured. As already discussed for the plasma fre- 
quency measurement, we partitioned the data — eigen- 
mode resonances are plotted with full symbols, whereas 
open symbols correspond to parasitic resonances. Beside 
a certain ambiguity of the data selection process, a good 
agreement with theory can clearly be seen. 

The measurement results of the more interesting AFM 
case are shown in Fig.^c). Here, injector currents were 
reset after each switching current measurement to ensure 
a well defined initial molecule state. Also, only eigen- 
mode resonances are plotted. Again, there is a good 
agreement between experiment and theory. It is inter- 
esting to note, that for k > tt almost no signature of 
the direct vortex state (-K, k) can be seen. Instead, the 
measured resonances correspond to the complementary 
state (27r — k, k — 2n), which has a higher eigenfrequency 
wq- As mentioned in Sec. Ill Bl for ir < k < k c both the 
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FIG. 7: Comparison of extrapolated zero bias eigenfrequency 
(symbols) with the numerically calculated ujq (k, 7 = 0) de- 
pendence (lines) of a vortex molecule in sample 1. (a) FM 
order, (b) AFM order. Dashed lines correspond to ideal dis- 
continuities, whereas solid lines take the finite injector size 
into account. 



direct and the complementary molecule state are stable, 
yet have different energies and different critical currents, 
cf. Fig. [5jd). In case of sample 1, k c ~ 1.37T. When a 
bias current is applied that exceeds 7 C of the direct state, 
the molecule flips into the complementary state by ex- 
change of one fluxon between the vortices. In contrast 
to the FM case however, here this transition is not auto- 
matically followed by the switching of the junction into a 
non-zero voltage state. So, as the experimental data in- 
dicates, already for k = I. In the vortex flipping usually 
remains undetected and only the more stable complemen- 
tary state is measured. 

Using Eq. ([3]) we fitted the measured resonance peaks 
(full symbols) in Fig. [5] leaving 7 c (k) and u}q(k) at zero 
bias as free parameters. Figure [7] shows the obtained 
values of ljo{ k ) and compares them with theory. As ex- 
pected from Fig. [6] we find very good agreement. As 
the finite injector size influences the eigenfrequency spec- 
trum, Fig. [JJ also shows the theoretical dependences of 
ujo(k) for ideal discontinuities. Clearly, only for large dis- 
continuities, i.e. k > 7r, differences become significant. 

In sample 1 only the lowest eigenmode could be in- 
vestigated due to limitations of our rf-generator. To ex- 
perimentally investigate higher eigenmodes we measured 
sample 3, which is the longest J J and also has the low- 
est plasma frequency. Figures |5[a) and (b) show the ob- 
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FIG. 8: (Color online) Comparison of measured resonance 
peak positions (symbols) with the numerically calculated 
U)±(k, 7) dependences (solid lines). Data is shown for (a) FM 
and (b) AFM molecule configurations in sample 3. Same col- 
ors correspond to the same k value. 



FIG. 9: (Color online) Comparison of the resonance peak 
positions of the lowest mode for annular LJJs of different 
lengths with two fractional Josephson vortices: experiment 
(open symbols), numerical simulations (solid lines), (a) FM 
(-0.6tt, -0.67r)-molecule. (b) AFM (-0.6tt, +0.6tt)- molecule. 



served resonance peak positions of the measured FM and 
AFM configurations. Parasitic resonances are not plot- 
ted. Interestingly, in the FM configuration both eigen- 
modes u>± can clearly be seen, yet only the lowest mode 
is present in the AFM state. Further measurements con- 
firmed, that for all symmetric FM states we selected, the 
observed resonances clearly showed two modes coinciding 
with the theoretical dependence of ui±, yet only a single 
mode for the AFM states. Even more surprising, mea- 
surements on two other samples of proper length (not 
shown) revealed no trace of u>\, neither in AFM nor the 
FM state, indicating some general problem with the exci- 
tation of higher modes — sample 3 being the exception. 
So far it is unclear to us, whether this behavior is owed 
to the particular way we apply microwaves to the JJ, 
i.e. by capacitivc coupling to the bias line, or to a more 
fundamental reason. 

As has been shown in Fig. [3] coupling between the 
vortices in a molecule can be adjusted by changing the 
distance a between them — the smaller is the distance, 
the larger is the coupling (frequency splitting). Clearly, 
compared to the FM configuration, the deviation of the 
frequency of the lowest mode from the single vortex 
frequency @ is much more pronounced in the AFM 
case. The configuration related, different dependence of 
ujq on a can also be seen in Fig. |Hl where we directly 



compare ^0(7) of the FM and AFM ordered molecules 
(pi, pi) = (-0.67T, -0.6tt) and (pi, p 2 ) = (-0.6tt, +0.67t) 
for three different vortex distances a. As expected, the 
eigenfrequency luq of the AFM state is much more sensi- 
tive to a than of the FM state. Finally, we used Eq. ([3]) 
to obtain Wo,(i)(^) at 7 = from measurements, which 
are shown in Fig.H 



IV. CONCLUSIONS 

We experimentally investigated the interaction of frac- 
tional Josephson vortices in a symmetric two-vortex 
molecule by measuring the oscillatory eigenmodes of the 
molecule spectroscopically. We observe a splitting of the 
single vortex eigenmode due to vortex-vortex coupling. 
The emerging molecule modes correspond to an in-phase 
and an out-of-phase oscillation of the individual vortices. 
The measured dependence of the eigenfrequencies of the 
molecule on applied bias current, topological charge p of 
the vortices and vortex ordering are in good agreement 
with numerical calculations. We find that with decreas- 
ing vortex distance coupling and frequency splitting in- 
creases. Depending on the configuration of the molecule, 
i.e. parallel or antiparallcl vortex order, spectra show 
significantly different behavior, which can, in principle, 
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